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Abstract. By subtracting from a recent lattice measurement of the thermal vector-current correlator the 
known 5-loop vacuum contribution, we demonstrate that the remainder is small and shows no visible short- 
distance divergence. It can therefore in principle be subjected to model-independent analytic continuation. 
Testing a particular implementation, we obtain estimates for the flavour-diffusion coefficient (2nTD > 0.8) 
and electrical conductivity which are significantly smaller than previous results. Although systematic errors 
remain beyond control at present, some aspects of our approach could be of a wider applicability. 



1 Introduction 

If a plasma is subjected to a perturbation whose wave- 
length is much longer than the typical mean-free path, 
then decohering scatterings take abundantly place within 
that perturbation, and its macroscopic physics should 
be classical in nature. For instance, if we imagine a 
hadronic "jet" with a width of a fermi or more, placed 
within a strongly interacting medium at a temperature 
T 3> 200 MeV, in which the characteristic distance scale 
is l/(nT) <C fm, then we expect to be able to de- 
scribe the gross features of the jet's subsequent behaviour 
within classical hydrodynamics. The main role of quan- 
tum physics is that the classical description involves pa- 
rameters, called transport coefficients, which need to be 
matched to the fundamental theory, in order to correctly 
capture the dynamics. 

As a particular case, let us consider a perturbation car- 
rying a specific net quark flavour. In the absence of weak 
interactions, there is a conserved current related to each 
flavour: d^Jf = 0, f = 1, Nf . The sum of all flavour cur- 
rents (divided by the number of colours, N c ) defines the 
baryon current, whereas a particular linear combination, 
weighted by the electric charges of all flavours, defines the 
electromagnetic current (denoted by J^ m ). 

Now, to define the flavour diffusion coefficient, denoted 
by D { , requires a specification of the classical description 
onto which to match. Let us assume that the perturba- 
tion is broad compared with 1/(ttT), and express Jf in a 
gradient expansion in ~ d^/irT. Apart from <9 M the other 
Lorcntz vector available is the four-velocity defining the 



fluid rest frame (m p 



1). The free coefficients al- 



lowed by Lorentz symmetry are the transport coefficients. 
In particular, we can expand 



where n f is the number density, and the transverse deriva- 
tive has been defined as 

= {rf v - u^ l u v )d v , rf = diag (+ ) . (1.2) 

The so-called Landau-Lifshitz convention is implied here, 
whereby in the rest frame of the fluid the zero component 
of Jf is nothing but the number density: u^Jjf = n f . 

To match for D { , it is helpful to specialize to simple 
kinematics. In particular, going to the fluid rest frame and 
imposing current conservation on eq. (1.1) yields 



d t m = D f V 2 n[ + C(V 3 ) 



(1.3) 



In this frame it is a textbook exercise to derive an expres- 
sion for Df in terms of 2-point correlators, which can then 
be matched onto quantum-mechanical expectation values; 
the result reads (cf. e.g. ref. [1]) 



1 

3x f 



lim > 



P»(w,0) 



(1.4) 



where Pii (u,0) = / tx e^{\ [J}(t, x), J}(0, 0)]) T is the 
spectral function related to the operators J}, and 



Xt = j|Vjf(i?(T,x)#(0,0) 



= /3 /(j°(0,x)J°(0,0) 



J' 1 



+ Dfd^n f + 0{d 2 ) (f fixed), (1.1) 



0=?, (1-5) 



is the "susceptibility" related to the conserved charge. It 
is important to stress that even though this way of deter- 
mining D f makes use of the fluid rest frame (reflected by 
the vanishing spatial momentum in eq. (1.4); in the follow- 
ing this redundant argument is suppressed) , the coefficient 
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itself is defined also for relativistic flow; the correspond- 
ing covariant form of the diffusion equation is eq. (1.1) 
together with d^Jf = 0. 

In the case of the electrical conductivity, a, a possi- 
ble definition is (J om ) = ctE, where E denotes an ex- 
ternal electric field. Let us denote D { = D fs + D fns , 
Xf = Xf. s + Xf, nS ' where D f , Xf „ are the contributions 
of the "singlet" or "disconnected" quark contractions, and 
D f ns , Xf ns those of the "non-singlet" or "connected" ones. 
One then obtains 



e 2 Xf 



[(E^f) Df, s + (Y,$) D f,™ 



f=i 



(1.6) 



where Q f denotes the electric charge of flavour f in units 
of the elementary charge e. 

Although a lattice Monte Carlo determination of D { 
and a is numerically very demanding [2] , a number of at- 
tempts have been launched in recent years [3-5] . In partic- 
ular, in ref. [5] a continuum extrapolation of the relevant 
(connected) Euclidean correlator was carried out for the 
first time, with a philosophy that analytic continuation 
from the Euclidean correlator to the Minkowskian spec- 
tral function should only be attempted on the continuum- 
extrapolated result. Subsequently various models were 
employed for the analytic continuation, based on a few- 
parameter ansatz for p u . In the present work we make use 
of the continuum-extrapolated Euclidean data of ref. [5], 
but analyze it in a different way, avoiding fits. 

The basic philosophy of our approach comes from 
ref. [6] , whose results were transcribed into a practical al- 
gorithm in ref. [7] . The main point is that in order to allow 
for an analytic continuation in principle, short-distance di- 
vergences need to be subtracted from the Euclidean corre- 
lator, such that Fourier coefficients can be determined (in 
fact the function should even be continuous [6]). An im- 
portant further insight comes from ref. [8], which showed 
that the ultraviolet (UV) asymptotics of the thermal con- 
tribution to the spectral function p u is such that it does 
lead to a continuous Euclidean correlator. Therefore, only 
the contribution of the vacuum p u needs to be subtracted. 
The final ingredient is that the vacuum p u can be ex- 
tracted from a recent 5-loop computation of the vector 
current correlator [9]. Implementing all these ingredients, 
we find surprisingly stable results which can be compared 
with other approaches, in order to obtain a rough impres- 
sion on the systematic uncertainties involved. 



2 Detailed setup 

We have in mind QCD with three massless valence flavours 
(N v = 3). The gauge field configurations were generated 
within pure SU(3) gauge theory in ref. [5], so the number 
of dynamical quarks is zero (Nf = 0). Moreover in ref. [5] 
only the "connected" or "non-singlet" contractions were 
evaluated. For N v = 3 this implies that we are technically 
considering electric charge diffusion and susceptibility; fol- 
lowing ref. [5] the corresponding coefficients are denoted 



by D = D f ns and x q = Xf ns - Otherwise we keep the no- 
tation of sec. 1 in the following, in particular continuing to 
employ Minkowskian conventions for the Dirac matrices. 
(Note that the physics of heavy flavour diffusion, referring 
to quarks with a mass M >• ttT, is quite different from 
that of light quark diffusion [10,11], and the two cases 
should not be confused with each other, even though a 
single notation D is often used for the diffusion coefficient; 
in particular, in the heavy quark case the extraction of D 
might be somewhat more robust than here [12,13].) 

To specify the observables, it is helpful to start with 
the Lorentz covariant form of the vector current correlator, 



3 „ 

Gy(r) =~^2 (OMOfox) WW)(0,0) 

,,— n -'x 



= -G 0Q (t) + G u (t), 



(2.1) 



where in the spatial part a sum over the indices is im- 
plied. Because of the projection to zero spatial momen- 
tum, G 0() (t) is actually r- independent, like in eq. (1.5); 
the value is denoted by G 00 (r) = x q T- It turns out that 
in the free limit, the spatial G u (t) contains the same con- 
stant, which then cancels in the sum of eq. (2.1) [14]. We 
denote this by Xq OC = N c T 2 /3 [15]. In the interacting 
theory, G 00 remains constant whereas G u gets essentially 
modified. 

The spectral functions corresponding to G v (t) and 
G u (t) are denoted by p v {oo) and P u {oj), respectively. If 
the spectral function is known, the Euclidean correlator 
can be obtained from 



raw 

Jo n 



cosh 



(M 



sinh 1 



(2.2) 



The basic issue is to what extent the inverse is true, i.e. 
information about the spectral function, particularly con- 
cerning the transport coefficient lim w _,. + p u (w)/uj rele- 
vant for eq. (1.4), can be extracted from a measured G u . 
An extensive review on the problems encountered and the 
methods currently available can be found in ref. [2]. 

We note in passing that the relation in eq. (2.1) im- 
plies a corresponding relation of the spectral functions, 
p v (cu) = Pa(w) — nx q ujS(uj). However the transport coef- 
ficients extracted from p v and p u as lim w _ >0 + /°( w )/ w are 
identical. In fact in the rigorous algorithm of ref. [6] the 
r- independent mode gets explicitly projected out. 

Now, in ref. [8], the UV asymptotics of thermal spec- 
tral functions were analyzed with Operator Product Ex- 
pansion methods. In particular, it was shown that ther- 
mal corrections to p u (tv) decrease at large frequencies as 
~ T 4 /u 2 . When taken together with eq. (2.2) this state- 
ment, valid beyond perturbation theory, implies that the 
thermal part of p h {oj) yields a contribution to G u {t) 
which is integrable even at r = 0, i.e. remains finite at 
short distances. In contrast, the vacuum spectral function 
yields a contribution diverging as ~ 1/r 3 . In order for the 
rigorous analytic continuation of ref. [6] to be applicable, 
the divergent part needs to be subtracted [7]. 
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In the literature, different strategies have been pursued 
for the subtraction. In particular, in refs. [16,17], analo- 
gous (but different) correlators were considered, and the 
idea was to compute the whole G(r), including thermal 
corrections, up to next-to-leading order (NLO), i.e. 2-loop 
level, or 0(a s ). Although the same strategy would also be 
possible for the vector correlator, by making use of clas- 
sic results for the NLO thermal spectral function [18-20], 
our goal here is to probe a different strategy. Namely, the 
vacuum subtraction is handled with much higher preci- 
sion than NLO, by making use of the fact that in vacuum, 
Pa(ui) is known up to 5-loop level, or 0{a A s ) [9]. In con- 
trast, the thermal part is handled with lower precision, 
only at leading order, leaving all other thermal effects to 
be taken care of by the non-perturbative numerical treat- 
ment of the remainder. 

Before proceeding it is important to underline once 
more the implications of the asymptotic behaviour ~ 
T 4 /uj 2 [8]. In particular, the Lorentzian form 



p£ L) H S 3D Xq 



tor) 



UJ 2 + I] 2 



(2.3) 



is sometimes used for modelling the transport peak. How- 
ever, this shape can be correct only at small frequen- 
cies [21]; at large frequencies to 77 it decays as ~ 
3Dx q V 2 1 ' w 7 which is slower than the mentioned ~ T 4 /uj 2 
(and does not allow for a negative sign which is also a 
possibility [8]). Most significantly, the Euclidean correla- 
tor obtained by inserting p^) m ^ G e q (2.2) diverges at 
t <C p. So, if used as a part of a fit ansatz for all uj, 
this function may pick up an incorrect overlap on vacuum 
contributions, which could lead to an overestimate of D. 



3 Vacuum spectral function 

We now turn to ref. [9] and specify the 5-loop vacuum 
spectral function. Following the conventions of ref. [22], 
the coefficients of the /3-function are defined according to 



d t a s = -(p a 2 + /3 ia 3 s + p 2 a\ + ...), (3.1) 



where 



a s = , t = ln( -j3- ] , (3.2) 



and, for iV c = 3 (cf. ref. [23] ; we only need terms up to the 
3- loop level here), 



11 JV f 

130 T ~ T ' 

. 51 19iVf 
Pi — - 



02 



8 24 ' 
2857 5033-/V f 325iV f 2 
128 1152 + 3456 



(3.3) 
(3.4) 
(3.5) 



The scale parameter in eq. (3.2) represents an in- 
tegration constant and, as usual, is chosen so that the 



asymptotic (r> 1) behaviour reads 

1 ftlni , p 2 (ln 2 t-lnt-l) + p 2 p /l 



Pot Pit 2 



■ + ■ 



m 3 



(3.6) 

Then, the main result can be obtained from the coeffi- 
cients Ilf, i = 0, 4, in eqs. (6)-(9) of ref. [22]. In vacuum, 
the quantity that we are interested in can be written as 



% > 2 



(3.7) 



where 



R{uj 2 ) =47rImfi7 + 77 1 a s + 77 2 a 2 + ..\ ; (3.8) 

L ~ J L -> i+i-!T 



and 



Writing 



(3.9) 



'3,1 

r 3,2 
r 4,0 



R{u?) = r + r lfi a s + (r 20 + r 2l i) a 2 s 

+ (r-4,o + r 4 ,i i + r ia I 2 + r 4 , 3 £ 3 ) a A s + 0(a 5 8 ) , 

(3.10) 

the generalization of eq. (5) of ref. [9] to I ^ reads 

r = 1.0000000 , 
r 10 = 1.0000000 , 
r 20 = 1.9857074 - 0.1152954^ , 
r 2 1 = 2.7500000 - 0.1666667iV f , 
r 3 Q = -6.6369356 - 1.2001341^ - 0.00517836Af , 
r, , = 17.296391 - 2.0876938iV f + 0.03843180Af , 
7.5625000 - 0.9166667iV f + 0.02777778Af , 
-156.60811 + 18.774765ATf - 0.79743434^ 

+ 0.021516105^ f 3 , 
r iA = -7.1166367 - 15.565615AT f + 0.83393599^ 

+ 0.002589178AT f 3 , 
r 4 2 = 88.878862 - 16.175418^ + 0.81239907^ 

- 0.009607950AT f 3 , 
r 4 3 = 20.796875 - 3.7812500AT f + 0.22916667Af 

- 0.004629630iV 4 3 . 

(3.11) 

As already mentioned, the lattice simulations of ref. [5] 
were for quenched QCD (iVf = 0) and only evaluated the 
"connected" quark contraction. In the language of refs. [9, 
22] the latter corresponds to "non-singlet" contributions, 
which are the only ones included in the results above. 
Therefore, eqs. (3.10), (3.11) with N { = can directly 
be used for the analysis of the lattice data of ref. [5] . (Ul- 
timately singlet contributions will need to be included as 
well, and on the perturbative side progress in this direc- 
tion is being made, cf. e.g. ref. [24].) 
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Fig. 1. Five orders of Gi 4 (r), from eqs. (2.2), (3.10), (4.2), nor- 
malized to the free result from eq. (4.3). (The gauge coupling 
has always been evaluated with 3-loop running.) 



4 Imaginary-time correlator 

Given the spectral function from eq. (3.7), re-interpreted 
as p Ul we could insert it into eq. (2.2) and compute the cor- 
responding imaginary-time correlator. There is the prob- 
lem, though, that technically arbitrarily small values of to 
contribute to the integral, but for those the determination 
of a s becomes ill-defined. However, any modification of p u 
in a finite range < uj < cj max yields a contribution to 
Gjj(r) which does not diverge at small t and can thus be 
taken care of together with the rest of the remainder. So, 
we cut the smallest frequencies off from a s by defining 



/x re f = max(7rT, u>) 



(4-1) 



and by choosing p, = /l rc f in the evaluation of G u (t). In 
addition we modify eq. (3.7) by accounting for the leading- 
order thermal corrections: 



ill 



l-2n F 



(I)] i^ + TTX^M , (4.2) 



where n F is the Fermi distribution. Following ref. [5] the 
temperature is set to T = 1.45T C , and following e.g. 
ref. [25] we take T c = 1.25/k^s (uncertainties related to 
this choice are discussed below). Like in ref. [5] the results 
are normalized to the expression [26] 



duj ( 3w 
2^ 117 



2 . 



1 - 2n ¥ 



(1) 

cosh (j^ ~ T ) 



ttT 2 ui5(uj) 



sinh^ 



6T C 



tt(1 - 2tT) 



1 + cos 2 (2tttT) 
sin 3 (27rrT) 



2 cos(2tttT) 
sin 2 (27rTT) 



(4.3) 



N = 24, o=10 ,/ =10 

max. 



<Q- 




Fig. 3. A test of the method of ref. [6] with the example in- 
troduced in ref. [7], employing an improved prescription. Here 
p refers to normalization by appropriate powers of T; and "in- 
put" denotes the correct result, known in this case. The result 
should be compared with fig. 4(right) of ref. [7], where an iden- 
tical data set lead to considerably more spread. 



In fig. 1, five subsequent orders of the perturbative result 
are shown; in fig. 2 the 5-loop perturbative result is com- 
pared with the continuum-extrapolated data from ref. [5] . 



5 Analytic continuation 

In order to estimate D from eq. (1.4), data for the Eu- 
clidean G H (t) should be fed into an analytic continuation 
prescription, yielding the corresponding p u {uj). The func- 
tion G^ rt based on eq. (4.2) yields a vanishing D, because 

(T) 

Pli has a vanishing slope around the origin. Therefore we 
can equally well apply analytic continuation to the differ- 
ence 

G-f " Gl ert . In fact, not only are we allowed to 
do this, but we probably must do this, given that to our 
knowledge mathematically justified analytic continuation 
to Minkowskian signature has been worked out only for a 
function which has no singularity at small r [6]. 

After the subtraction of the singular terms, we make 
use of the method of ref. [6], as implemented in ref. [7]. 
In the meanwhile we have even found a significant im- 
provement over the original implementation. The general 
method is based on determining the Fourier coefficients 
from the Euclidean data G(r), which can be done nu- 
merically, and using these in order to construct an expan- 
sion of the real-time dependence of the correlator in terms 
of Laguerre polynomials, with coefficients denoted by at. 
The physical real-time correlator must vanish at infinite 
time separation [27], which corresponds to 5^=0* ae = 0- 
In practice, however, the correlator does not vanish ex- 
actly; choices of £ max for which it vanishes approximately 
were dubbed "windows of opportunity" in ref. [7] . Within 



Y. Burnier, M. Laine: Towards flavour diffusion coefficient and electrical conductivity 



o 



1.06 



1.04- 



1.02- 



1.00 




0.2 



0.0 - 



o 

I 



-0.2 



-0.4 



0.0 



T= 1.45 T ,T = 1.25 A™ 

c' c MS 








o.i 



0.2 









= 0.900 T 2 


o x q 


= 0.897 T 2 




= 0.894 T 2 




0.3 


0.4 



0.5 



XT 



Fig. 2. Left: Comparison of the O(a^) result from fig. 1 with lattice data from ref. [5], multiplied by Xq/T 2 = 0.897. In the 
0(<is) result, the renormalization scale has been varied within the range ft — (1...5) x ft TO f (grey band). Right: The difference of 
the lattice and perturbative results, with Xq varied within its uncertainties. Errors are shown for Xq/T 2 = 0.897 only. 



these windows, the small remaining asymptotic value was 
subtracted from the correlator by hand, in order to allow 
for a Fourier transform to Minkowskian frequency space. 
We have now replaced the subtraction by another proce- 
dure; namely, the coefficient a e for which the asymp- 
totic value first crosses zero, is redefined to be 



E 

£=0 



(It 



(5.1) 



It turns out that this way the dependence on £ max is much 
milder than with the original procedure; in fact one does 
not even need to be close to a "window" but a plateau in 

=o* a ? wou ld be sufficient as was envisaged in ref. [6] . 
The improvement is illustrated in fig. 3, obtained for the 
same data set (with the same simulated errors) as fig. 4 
of ref. [7]; the results are substantially more stable, and 
equally close to the correct value ("input"), underestimat- 
ing it by ~ 25%. (We have shown results for N = N T — 24 
to allow for a direct comparison with ref. [7]; the case 
N = 48 corresponding to the resolution of ref. [5] shifts 
the results by 2-3% in the correct direction.) 

Encouraged by this success, we have applied the al- 
gorithm to the difference G^ lff = G-f u — G^ Glt shown in 
fig. 2 (right) (the corresponding spectral function is de- 
noted by pfi 3 )- 1 At very short distances no lattice data 
exists, so the difference has been kept fixed at its value at a 
chosen r = T m j n (note that there is no reason for the differ- 
ence to vanish). Various sources of systematic errors have 



1 The data for G^ lff as well as for a corresponding G v *,<,,, 
be obtained from the authors on request. The two differ by a 
constant mode; we stress that in the algorithm of ref. [6] or in 
the fit described in footnote 2 an exactly constant mode has 
no effect on the spectral function. 



-.diff 



been probed. We have checked that variations of the renor- 
malization scale, as indicated in fig. 2(left), and variations 
of T c /A^ within a range 1.25 ± 0.10, consistent with e.g. 
ref. [25] (but also with earlier works, cf. sec. 4.2 of ref. [28]), 
have an effect smaller than variations of X q /T 2 = 0.897(3), 
to which the continuum-extrapolated results were normal- 
ized in ref. [5]. The errors from the latter variation are 
shown in fig. 4(left). As can be anticipated from fig. 1, 
using 3 or 4-loop vacuum results would only lead to minor 
changes. Another source of systematic errors is the contin- 
uum extrapolation of ref. [5] ; it might be prudent not to 
make use of results below tT ~ 0.20 [5]. In fig. 4(right) we 
show the corresponding effects; in the following we restrict 
to r m i n T = 11/48 which lies within a stable range. 

It remains to consider statistical errors. The errors as 
shown in fig. 2 are strongly correlated, but only available 
for each tT separately. We have then shifted the whole 
function upwards or downwards by the errors. It is im- 
portant to stress once again that the algorithm of ref. [6] 
exactly projects out any constant contribution (i.e. Mat- 
subara zero mode), so a uniform shift has no effect; the 
difference comes from the change in the shape. Results 
are shown in fig. 5; the variation from this rough imple- 
mentation is smaller than systematic errors related to the 
analytic continuation. 



6 Conclusions 

Based on figs. 4, 5, and folding in an estimate of a down- 
ward systematic error as suggested by fig. 3, we estimate 

~ diff (w) 



lim 



lim 



Pi 



> 0.3 T 



(6.1) 



Taking into account a factor 2 difference in the nor- 
malization of p u , the corresponding estimate was cited 
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Fig. 4. Left: Dependence of the spectral function on Xq/T 2 , for a fixed r m - m T — 11/48. Right: Dependence of the spectral 
function on r mul T, for a fixed Xq/T 2 = 0.897. 
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Fig. 5. Results obtained by modifying lattice data by adding 
or subtracting its statistical errors. 

as (1...3)T in ref. [5], i.e. 3-9 times larger than our 
lower bound. (In ref. [29] only the lower limit was quoted, 
Dxq — 0.33T. These estimates rely on eq. (2.3) in com- 
bination with non-logarithmic modellings of the UV part 
of p u 2 ) From eq. (1.4), inserting Xq = 0.897T 2 , eq. (6.1) 

2 If the remainder from fig. 2(right) at tT > 11/48 is sub- 
jected to a 3-parameter fit to eq. (2.3) plus a r-independent 
constant, and errors are treated as uncorrelated, we obtain 
D Xq ^ 0.30T, 7j ~ LIT, with xVd.o.f. ~ 0.006, numbers 
comparable with ref. [29]. The corresponding transport peak 
is markedly narrower than those in figs. 4, 5. Yet, as men- 
tioned, the associated Gf? S diverges at small r, whereas ours 
stays finite, as required by ref. [8] . If we cut off the large fre- 



yields 

2nTD > 0.8 . (6.2) 

The corresponding electrical conductivity from eq. (1.6) 
evaluates to a > 0.07 e 2 T for iV v = 3, N f = 0. 

These numbers are substantially smaller than those 
of the leading-order weak-coupling expansion [30,31], but 
they are intriguingly close to the AdS/CFT suggestion 
2irTD = 1 [32]. However, the transport peak could be ex- 
tremely narrow [33] and therefore our results should be 
interpreted as lower bounds, as has already been indi- 
cated by the notation. Continuum results are needed in 
the full r-range, and the resolution needs to be gradually 
increased, both in terms of statistical precision as well as 
in terms of N T , in order to see whether the results stay 
put or show a slow evolution. (Unfortunately the func- 
tional dependences of the analytically continued results on 
statistical variance and N T are not easily extracted [6].) 

In any case, the principal feasibility of carrying out 
a model-independent short-distance vacuum subtraction 
before attempting any analytic continuation or spectral 
modelling, thereby removing harmful UV contamination 
from the signal (cf. fig. 2), has hopefully become clear. 
This general philosophy can perhaps be applied to other 
correlators as well, for instance those related to compo- 
nents of the energy-momentum tensor, even if in that case 
temperature-dependent subtractions are needed in addi- 
tion to the vacuum one [8] . 



quencies responsible for the divergence, for instance by defining 
(w) = p^ (w)/ cosh(2^rj0, then the fit result moves in our 
direction: D Xq ^ 0.16T, r\ ~ 2.0T, with xV d -°L ^ 0.005. To 
resolve the correct physics it hence appears important to reach 
a good resolution for the continuum extrapolation of Gfl S also 
at small r, verifying that it saturates to a constant value there. 
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